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We study the nonequilibrium interplay between disorder and interactions in a closed quantum 
system. We base our analysis on the notion of dynamical state-space localization, calculated via 
the Loschmidt echo. Although real-space and state-space localization are independent concepts 
in general, we show that both perspectives may be directly connected through a specific choice 
of initial states, namely, maximally localized states (ML-states). We show numerically that in the 
noninteracting case the average echo is found to be monotonically increasing with increasing disorder; 
these results are in agreement with an analytical evaluation in the single particle case in which the 
echo is found to be inversely proportional to the localization length. We also show that for interacting 
systems, the length scale under which equilibration may occur is upper bounded and such bound 
is smaller the greater the average echo of ML-states. When disorder and interactions, both being 
localization mechanisms, are simultaneously at play the echo features a non-monotonic behaviour 
indicating a non-trivial interplay of the two processes. This interplay induces delocalization of the 
dynamics which is accompanied by delocalization in real-space. This non-monotonic behaviour is 
also present in the effective integrability which we show by evaluating the gap statistics. 

PACS numbers: 72.15.Rn; 71.30.+h; 05.30.Fk 


I. INTRODUCTION 

The physics of transport in electronic systems is a cor¬ 
nerstone of our technologies and yet the general prob¬ 
lem of non-equilibrium transport seems far from being 
completely understood. Significant steps in this direction 
have been taken starting with Anderson localization (or 
absence of diffusion(P. It is well established that in the 
absence of interactions the energetic mismatch between 
neighboring sites in a lattice may completely prevent 
transport by localizing all single-particle states. In this 
scenario arbitrary weak disorder localizes one and two- 
dimensional systems, producing ideal insulators. Much 
less is known about the effects of the interplay between 
interaction and disorder. In fact, despite intensive studies 
over decade^^HH^ the influence of electron-electron inter¬ 
actions on the transport in disordered electronic systems 
is still a challenging problem. 

More recently it has been theoretically shown that the 
insulating phase survives the inclusion of a small amount 
of interaction in a disordered system although, in this 
case, one ca n induce transport by increasing the tem¬ 
perature 1^. This investigation has given rise to the 
emerging field of many-body localization (MBL), which is 
considered a dynamical transition; unlike an equilibrium 
(conventional) phase transition, in a MBL transition the 
properties of the entire set of eigenstates change and one 
cannot restrict the analysis to the ground-state physics. 
This dynamical localization has important consequences 
on time evolution, more specifically on the capacity of the 
system to self-equilibrate, that is, to act as a reservoir for 
its subsystems. Localized states may prevent equilibra¬ 
tion, since local perturbations do not diffuse throughout 


the system - such systems are called non-ergodi(P. 

In the last few decades we have witnessed both exper- 
imentaPtii^ and theoreticaP^^li^ advances in the physics 
of one-dimensional quantum systems. In particular, af¬ 
ter Basko’s seminal worlP many groups have focused 
on the non-equilibrium localization of interacting parti¬ 
cles. In these works the MBL transition has been stud¬ 
ied through different approaches, using, for instance, en¬ 
tanglement properties or eigenstate statistics. Numerical 
studies have shown that the spectral statistics of a one 
dimensional system crosses over from those of orthog¬ 
onal random matrices in the diffusive regime at weak 
disorder to Poisson statistics in the localized regime at 
strong disordeifi^. Concomitantly, dynamical features of 
disordered and interacting systems were observed focus¬ 
ing on multi-partite correlationJi^ and on the entangle¬ 
ment entropy. It was shown that the interplay between 
interactions and disorder has a strong influence on the 
long-time behaviour of the entropjff^*® later on it has 
been concluded that such signatures are just a character¬ 
istic trait of interacting systems, since these dynamical 
tendencies also exist in the absence of disorder and inho- 
mogeneitie^^. An analysis of the spatial behaviour of the 
entropy (von Neumann entropy as a function of partition 
size) should be more appropriate to address localization, 
since it saturates when the partition size reaches the lo¬ 
calization lengttP^. However this is a difficult approach, 
since the investigation has to be applied to large chains. 
In fact, the definite evidence of a MBL transi tion re mains 
elusive despite creative and innovative effortdue 
to technical difficulties, such as finite size limitation of 
numerical calculations. 

In this work we study the dynamical state-space local- 
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ization of an initial state that is maximally localized in 
real-space and evolves in time with an interacting and 
disordered Hamiltonian. In this way we are able to de¬ 
scribe the non-equilibrium localization effects of disorder 
and interaction in a well defined way which is accessible 
for small chains. Furthermore, we propose a direct re¬ 
lation between dynamical localization in the state-space 
and the more standard real-space localization perspec¬ 
tive. We show that when acting alone both disorder and 
interaction progressively and monotonically localize the 
system. However, when combined, they may lead to de- 
localization. 

The paper is organized as follows. Section |H] is dedi¬ 
cated to the model. We study a one-dimensional (1-D) 
spinless fermion system with nearest-neighbor Coulomb 
interaction and on-site disorder. The quantity we 
choose to address dynamical state-space localization (the 
Loschmidt echo) is described in section HI as well as the 


connection with real-space localization. We present our 
numerical results in section ITVl The last section is dedi¬ 
cated to our comments and conclusions. 


II. MODEL 


The disordered Anderson-Hubbard model is a standard 
Hamiltonian used to describe the competition between 
kinetic energy, electron-electron interaction, and disor¬ 
der. The fully polarized Anderson-Hubbard model, or a 
disordered spinless fermion model, can still describe this 
competition with less computational effort. In the latter, 
we only retain the nearest-neighbor Coulomb interaction 
since a site cannot be doubly occupied. The Hamiltonian 
is then given by 


^ = E 

(*d> 



+[/(n,-l/2)(n,-1/2)]. 


( 1 ) 


The first term describes the kinetic energy, with J being 
the tunneling between nearest neighbors; the second de¬ 
scribes the local energy of the ith site, with Ci taken from 
a uniform distribution [—IT, W]; finally, the last term rep¬ 
resents the nearest-neighbor Coulomb repulsion, with U 
being the interaction strength, q and are the destruc¬ 
tion and number operators acting on site i, respectively. 
Throughout this paper we use J = 1 as the unit of en¬ 
ergy. This model can be mapped to the XXZ spin-1/2 
chain through a Jordan-Wigner transformatiorPS. 

In this work we study small chains with N = 
10,12,14,16 sites and solve the Hamiltonian by exact 
diagonalization. Although the study of small chains is 
not ideal to describe phase transitions (defined in the 
thermodynamic limit), it can still produce relevant phys¬ 
ical insights about the interplay between interactions and 
disorder. Furthermore, small systems can be realized in 
chains of trapped ionP^ and cold atom^iS! a nd have been 
the subject of many recent theoretical studie^ ^T^b^ l ^^ l ^ ^. 


III. NONEQUILIBRIUM LOCALIZATION 

We choose to describe the competition between disor¬ 
der and interactions through a state-space localization 
perspective. The dynamical state-space localization can 
be described by the inverse participation ratio (IPR)P^ 
of the time average ensemble. The IPR is calculated 
through the purity of the time-average ensemble, which 
gives us an estimate of the inverse of the effective length 
or area covered by the dynamics: 

/Pi? = tr{p2}, p= lim ^ f \-^t){'i’t\dt, (2) 

T-!-oo 1 Jq 

where Idti) represents the dynamically evolving state at 
a time f, starting from an initial state l^/o). The IPR is 
bounded between unity and 1/d^ with d being the site 
dimension and d^ the dimension of the state-space. If 
IPR = 1 the dynamics are maximally localized mean¬ 
ing that an initial pure state is unaltered throughout the 
evolution; if IPR = 1/d^ the dynamics are maximally 
delocalized in the sense that the time-average ensemble is 
the identity ensemble in which all states contribute with 
the same probability. 

The purity of the average ensemble can also be ex¬ 
pressed by the time-average of the Loschmidt echo Lit) 

1 r 

IPR = L = lim — dt L{t) (3) 

T^OO T Jq 

where L{t) = |(4'o|e“®^*|'I'o)p. This relation becomes 
useful in the case where the full diagonalisation or long¬ 
time dynamics are unreachable (which is the case exper- 
imentalljPl) and we may use the dynamical behaviour of 
the echo as a figure of merit. Furthermore, the echo has 
proven to be an important tool in characterizing several 
aspects of the closed system dynamic^SHMl, 

Considering a non-degenerate Hamiltonian, the above 
relation corresponds to 

L = Y,\{mo)\\ ( 4 ) 

E 

where \E) represents the Hamiltonian eigenstates. In this 
way, the time averaged echo is related to the fidelity be¬ 
tween the initial state and the Hamiltonian eigenstates. 

In this work we calculate the dynamic state-space lo¬ 
calization through eq. Q. In summary, the dynamics 
of an initial state allowed to evolve in time with a given 
Hamiltonian is localized in the state-space when the par¬ 
ticipation ratio is close to unity (1/L ~ 1), while it is 
delocalized if 1/L « meaning that the state covers a 
large portion of the state-space through the time evolu¬ 
tion. 


A. Connection between state-space and real-space 
localization 

The equilibrium real-space localization, as introduced 
by AndersoiP, is a single-particle effect in which the 
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wavefunction exhibits exponentially decaying tails. In 
this way, each particle can only be found, with consider¬ 
able probability, inside a finite localization length ^ from 
a given lattice site. The definition of real-space localiza¬ 
tion in a many-body system is more subtil and in gen¬ 
eral we would refer to a correlation length rather than 
a localization length. Here we relate state-space local¬ 
ization with real-space localization by investigating the 
participation ratio of the dynamics of a special set of ini¬ 
tial states, namely half-filled random product states, i.e. 
states of the type |'Po) = and = N/2, 

where \ni) = {|0), |1)} are the eigenstates of the number 
operator of the “i-th” site, and N is the number of sites 
in the chain. Note that such states represent electrons 
which are maximally localized in real-space, with a van¬ 
ishing correlation length (^ —>■ 0), and we will name their 
set {|iil 5 _>.o)}.We perform an average of L over M initial 
product states and disordered Hamiltonian distributions 

M 

n E E r < 1, 

i E' E 

( 5 ) 

where (( )) denotes the total averaging, {e} represents 
the set of drawn on-site energies and = -^({c}) 

is the probability of a disorder configuration. 

In the extreme case in which a given product state 
allowed to evolve in time with a Hamilto¬ 
nian whose eigenstates {|Tl{«o)} are strongly localized 
in real-space (i.e. very weakly correlated), it will nec¬ 
essarily have a high projection on one of them, i.e. 

~ ^ee'- la this way, the state almost 
does not change in time, covering a small portion of the 
Hilbert space during the time evolution (state-space lo¬ 
calization), resulting in !/((£)) « I. On the other hand, 
if the eigenstates of H are {|T^{>i)} delocalized in real- 
space, |i?^_,,Q) will have a small projection all of them, 

l(^5>i|^.c_>o)P ^ 1- Therefore, the state covers a large 
portion of the Hilbert space during the dynamics, re¬ 
sulting in a large participation ratio 1/((L)) :g> 1 (dy¬ 
namical state-space delocalization). In general, we ex¬ 
pect that, when using random product states (which are 
maximally localized with vanishing correlation length) as 
initial states the dynamical localization becomes closely 
connected to the real-space localization of the Hamilto¬ 
nian eigenstates. Therefore, we expect a functional de¬ 
pendence {{L)) oc /(^), such that /(^) is a decaying func¬ 
tion of 

It is important to emphasize that the direct relation 
between real-space and state-space localization depends 
on the choice of initial state. Indeed, there are cases 
where the states are delocalized in real-space but local¬ 
ized in state-space. The simplest example is to consider 
initial states that are (at least approximately) eigenstates 
of a delocalized Hamiltonian. In this case the state would 
have an extremely long correlation length, however, the 
echo would remain (almost) unaltered throughout the 
evolution. 


B. Average participation ratio of maximally 
localized states and single particle localization 

One example where /(^) can be explicitly obtained 
is that of noninteracting particles in disordered systems 
at low density. In this regime we may approximate 
the problem by the single particle analysis which al¬ 
lows for a direct analytical treatment of the quantum 
states. The presence of scattering centers in a poten¬ 
tial induce bound eigenstates which are exponentially 
decaying around localization centers at X as |i?j(A)) = 
/ with a localization length ^ that is 

proportional to the depth of the potential. 

In our analysis we always initiate the system in a max¬ 
imally localized state which in the single particle case 
translates as \x') = lim^_,.o \E^{x')). Let us now evaluate 
the echo dependence on the localization length for this 
particular case: 

((!)) (X ly dx'l dX\{E^{X)\x')f = 1, (6) 

showing that the average echo is inversely proportional to 
the localization length and /(^) = | as we expected from 
the reasoning of the last section. Therefore, if the Hamil¬ 
tonian has very deep effects inducing a vanishing local¬ 
ization length the dynamics is localized in state-space. In 
the other limit when the Hamiltonian is almost clean with 
very large localization length the dynamics is highly delo¬ 
calized. Hence, we establish an explicit relation between 
the state-space dynamical localization and the real-space 
localization. Certainly, such simple functional relation is 
not necessarily expected for interacting particles. 

C. Average participation ratio and equilibration 

The previous reasoning can also be seen from the equi¬ 
libration perspective that allows for the analysis of a 
broader class of Hamiltonians which may account for 
interactions. The fact that the dynamics of the initial 
states (which are chosen to be localized in real-space) are 
localized in the Hilbert space implies in the breakdown 
of the ergodic hypothesis; the value of a time averaged 
observable depends on the initial state and the system 
does not equilibrate. The absence of equilibration has 
been used as an indicator of real-space many-body local¬ 
ization in both experimentaP and theoreticaP^ studies. 

More specifically, if one splits a system of N sites in two 
parts {S and R) and chooses two different initial product 
states |7rs)®|7rfl) and |7rs)®|7ri^) (with [tts) orthogonal to 
Itts)), equilibration guarantees that this choice becomes 
undetectable, for any practical purpose, if given enough 
time. This means, in particular, that the time average of 
the trace distance D{p,p) = — ps||i of the reduced 

evolved states of S, ps{t) and ps{t), should approach 
zero as t —)■ oo: lim2-_,,oo ^ D{psit), ps{t))dt — >■ 0. 
On the other hand, the distance remains large in the 
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absence of equilibration. In the appendix, we prove that 
if 1/ ((L)) < 1 + i5, with d <C 1 related to the size Ns of 
the subsystem S as then 

lim ^ f D{ps{t),ps{t))dt>l-6VS. 

T^oo T Jo 

It is thus shown that if the dynamics of the global sys¬ 
tem is localized in state-space then it is guaranteed that 
large enough subsystems composed of Ns sites do not 
equilibrate. This means that correlations do not extend 
over Ns sites and thus the correlation length should be 
much smaller than this subsystem. 


IV. NUMERICAL RESULTS 

We start by evaluating the participation ratio, 1/((L)), 
through eq. ([5) for the two simplest cases, that is, using 
a disordered Hamiltonian in the absence of interaction 
and also an interacting non-disordered Hamiltonian. We 
perform an average oven 10^ maximally localized initial 
states and, in the disordered case, for each initial state we 
randomly select each site energy with uniform probability 
between —W and W. In the absence of interactions the 
presence of disorder localizes the system progressively de¬ 
creasing its localization length. This result is represented 
in Figj^ by the monotonic decrease of the participation 
ratio as we increase disorder, which means the purity of 
the average ensemble only increases. In this case the de¬ 
crease of the participation ratio follows the reduction of 
transport due to Anderson localization. 

Similarly, for a non-disordered Hamiltonian, interac¬ 
tion progressively localizes the system. However the 
strongly interacting system is less localized than the 
strongly disordered system. The state-space localization 
analysis seems to indicate a partial transport suppres¬ 
sion due to interaction. This is in agreement with the 
equilibrium counterpart with a non-null conductivity in 
a system with large interaction but hnite temperatur^^. 

The difference between both interaction and disorder 
routes to localization shown in the inset of Fig. [^can be 
understood by applying standard time-independent per¬ 
turbation theory to the model Hamiltonian of equation 
Q. For a specific disorder realization, the disordered 
term of the Hamiltonian is diagonal in a basis of random 
product states and the spectrum of a finite system is 
non-degenerated. Thus, for strong disorder any product 
state tends to an eigenstate of the system. In the limit of 
non-interacting and strongly disordered one-dimensional 
systems all the eigenstates of the Hamiltonian Q are 
strongly localized; this result leads to 1/Z « 1. 

The electron-electron interaction term of the Hamilto¬ 
nian Q is diagonal in the same basis. Its eigenstates, 
however, are degenerated in sub-spaces of equal number 
m of nearest-neighbor filled sites, which we will refer to as 
Um subspaces. By including the hopping via degenerated 
perturbation theory we break the energy degeneracy of 
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FIG. 1: (Color online) Participation ratio (1/{(L))) of maxi¬ 
mally localized states evolving with a non-disordered Hamil¬ 
tonian as a function of interaction U (black circles) and with a 
non-interacting Hamiltonian as a function of disorder strength 
W (red squares). The inset shows a zoom for large interac¬ 
tion/disorder. Results obtained for a system of A = 12 sites 
through an average over 10'* initial states and disorder real¬ 
izations. 


these eigenstates and the new eigenstates of the system 
become superpositions of localized states belonging to the 
same Um subspace. The dynamics of a given initial prod¬ 
uct state with m nearest-neighbor filled sites is restricted 
to the corresponding Um subspace being, therefore, par¬ 
tially localized. That is why, in the strongly interacting 
limit, we still observe some residual dynamics. 

The competition between disorder and interactions 
leads to a challenging scenario which can not be de¬ 
scribed through perturbation theory, but the effects can 
still be observed in the participation ratio, as can be 
seen in Fig. Even though their individual natural 
effect is to localize the system, the combination of dis¬ 
order and interactions leads to nontrivial behaviour. In 
fact, the introduction of disorder in an interacting sys¬ 
tem (top panel) or of interactions in a disordered system 
(bottom panel) leads to delocalization. Such interplay 
continues up until the point at which one of the effects 
starts to dominate, i.e. when either disorder or inter¬ 
action becomes the dominant effect. More specifically, 
the interplay between interaction and disorder produces 
a tendency to delocalization when both U and W are 
of the same order as the hopping amplitude J. These 
results are imprinted in the non-monotonic behaviour 
of the participation ratio. In the thermodynamic limit 
this behaviour could lead to a reentrant phase diagram 
similarly to the infinite temperature equilibrium diagram 
outlined in|TT|and the non-disordered non-equilibrium di¬ 
agram outlined in 1421 noting that in the present case it is 
a genuine non-equilibrium effect of the interplay between 
interactions and disorder. 

The nonequilibrium interplay between disorder and in¬ 
teractions becomes more pronounced for larger systems. 
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FIG. 2: (Color online) Participation ratio, 1/{{L)), for (top) 
disorder strengths IF as a function of interaction U, and (bot¬ 
tom) different interactions [/ as a function of disorder W. The 
results were obtained averaging over different initial ran¬ 
dom product states and disorder realizations for a chain of 
= 12 sites. 


This is evident in Fig. where we show the rescaled in¬ 
verse echo for different system sizes as a function of the 
interaction strength. We scale the echo by the value it as¬ 
sumes when the interaction vanishes. In that case when 
{{L{U = 0)))/((L)) > 1 the interplay induces delocaliza¬ 
tion. In the top panel of Fig. we see that the region 
and magnitude above unity increase with the system size. 
From the bottom panel we confirm that this effect be¬ 
comes stronger for larger systems. We point out that 
the total mutual entropy has also been shown to have a 
non-monotonic behaviour^. This, once again, hints on 
the fact that the interplay is also manifested in the spa- 
cial correlations and that, indeed, defined with respect to 
maximally localized states, the average echo is a useful 
quantifier of localization out of equilibrium. We note that 
all the curves for different N cross at the same point and 
we could argue that such crossing could indicate a phase 
transition in the same spirit that has been previously 
done in the literature for other quantities. However, we 
feel that such a claim is unsubstantiated in the sense that 
such chain sizes are too small for extracting information 



FIG. 3: Scaled inverse echo, (((Z(t/)))/((Z(t/ = 0)))) ^ for 
W = 0.5 (Top) as a function of interaction strength and (Bot¬ 
tom) as a function of the number of sites. Results averaged 
over 5000 samples for N = 10, 2000 samples for N = 12, 1000 
samples for N = 14, and 500 samples for N = 16. 



FIG. 4: Inverse echo for W — 0.5 as a function of interac¬ 
tion strength for different magnetization subspaces. Results 
averaged over 1000 samples for N = 14. 


about the correlation length in the thermodynamic limit. 

Figures to exhibit results for a subspace of half- 
filled chains (zero magnetization subspace, in the spin 
scenario); in Fig. we show the inverse echo for dif¬ 
ferent filling sectors, analyzing how the results depend 
on the filling fraction, that is, the number of particles 
over the number of sites. Add one and two particles to 
a system at half-filling (corresponding to subspaces of 
magnetization m = 2 and m = 4, in spin units, respec- 
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FIG. 5: (Color online) Participation ratio for W = 1 obtained 
starting from specific half-filled initial product states. The 
solid black line corresponds to the data of Fig. i i.e., m 
is obtained through an average over 10^ initial product states 
and simultaneous disorder realizations. In the other curves, 
we start from the initial states specified in the legend and 
average L over 10^ disorder realizations. Results for N = 12. 


tively) progressively localizes the system, as indicated by 
the reduction of the participation ratio. The non-trivial 
interplay between interaction and disorder, which leads 
to a non-monotonic behavior of the inverse echo, is still 
present in the sectors we have analyzed, however it be¬ 
comes less pronounced. In fact, in the limiting case of 
fully occupied or fully empty system (maximal magneti¬ 
zation) disorder and interactions would not play any role 
in the participation ratio. 

Our results should be distinguished from the results 
in|43l in which a local echo is defined based on the magne¬ 
tization. Firstly the non-monotonic behaviour in|33|is al¬ 
ready present both at zero disorder (and varying interac¬ 
tion) or zero interaction (and varying disorder); secondly 
even though local probes are highly desirable experimen¬ 
tally, this one in particular is not a direct measure of lo¬ 
calization, although a reasonable indicator in some cases. 
We also point out that our findings for the Loschmidt 
echo differ significantly from those in|44|for several rea¬ 
sons: there, the authors use a class of initial states that 
are local representatives of an energy shell, while our 
choice is not based on energy but rather randomly se¬ 
lected from the entire set of non-entangled states. Their 
focus is also on the Heisenberg point while we sweep the 
interaction strength. 

As described above, we perform an average over dif¬ 
ferent initial random product states of full or empty 
sites, resulting in a half-filled chain (we remember that 
the model can be mapped onto a spin-half XXZ chain, 
and in the spin scenario this set of states is equivalent 
to random product states of spins up and down, with 
null magnetization). This set is a good representative of 
real-space localized states in general, which are strongly 
non-entangled. Our choice of initial states is essential 
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FIG. 6: (Golor online) Typical gap, given by the position 
of the maximum of the gap distribution, as a function of U 
(top) and W (bottom). The inset in the top panel shows 
the gap distributions for a system with W — 1 and different 
values of U, while the bottom inset shows the distributions for 
U = 1 and different disorder strengths. Results for A = 12 
and averaged over 5000 (top) and 20000 (bottom) disorder 
realizations. 


to the observation of our results and to the assurance of 
our interpretations, as confirmed by the results of Fig. 

Performing the same analysis described above, but 
starting with a specific local state, one could obtain dif¬ 
ferent results. Using a charge density wave state (CDW 
state, 1101010...)) as the initial state, for instance, we ob¬ 
serve a monotonic decrease of the participation ratio as 
we increase interaction in a disordered system (see Fig. 
[^. It happens because the CDW state, which is the 
ground state of the strongly interacting clean system, is 
not a typical localized state. It is contained in the Uq (no 
nearest-neighbor-filled sites) subspace, spanned by only 
2 of the 924 possible product states in a half-filled chain. 
Fig. 0 also shows that starting with a state that is a 
simple site-permutation of the CDW state is enough to 
recover the non-monotonic behaviour. 

To complete the analysis, we show that the interplay 
between disorder and interactions is also manifested in 
the gap distribution of the Hamiltonian, T’(A), where 
A = Ei+i — Ei and {Ei} are the Hamiltonian eigenval- 
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ues in ascending order. The gap distribution is commonly 
used as an indicator of integrability or chaos. Generi- 
cally, integrable system exhibits Poissonian gap statistics 
while chaotic systems exhibit Wigner-Dyson statistic^^. 
Furthermore, we expect localized regimes to resemble in¬ 
tegrable system and their statistics (since they have lo¬ 
cal constants of motion) while delocalized crimes are 
expected to resemble non-integrable systemPS A sub¬ 
stantial difference between these two regimes is that the 
maximum of the distribution moves closer to zero in case 
of localization and to positive values in case of delocal¬ 
ization. In Fig. I^we keep track of the typical gap, Aj^p, 
given by the position of the maximum of the distribution 
(examples of P(A) for disordered and interacting systems 
are given in the insets). This figure clearly shows that 
the gap statistics of a clean interacting system and of 
a disordered non-interacting system is always a Poisson- 
like distribution since they are integrable. On the other 
hand, the inclusion of disorder in an interacting system 
or the inclusion of interaction in a disordered system, 
leads to a non-monotonic behaviour of the typical gap, 
indicating that the nontrivial interplay between disorder 
and interactions drives the system to a more delocalized 
regime when W ~ U k, J. Atyp moves closer to zero as 
we increase disorder or interaction even further, recover¬ 
ing the localized regime. This analysis agrees with our 
previous results, that is, the maximum of Atyp as a func¬ 
tion of disorder or interaction is established for the same 
range of parameters in which we observe the maximum 
of the participation ratio (Fig. [^. 

V. CONCLUSION 

We have used the time-averaged Loschmidt echo of 
maximally localized states to characterize state-space lo¬ 


calization. We applied this concept to show that for a sin¬ 
gle particle the average echo is inversely proportional to 
the localization length, highlighting an intimate relation¬ 
ship between dynamical state-space localization and real- 
space localization. We have then extended the idea to 
many-body localization where we have derived a bound 
for subsystem equilibration based on this echo. Using 
this approach, we have shown that in a system featuring 
disorder and repulsive interaction, both localizing mech¬ 
anisms, the interplay between them can actually lead to 
delocalization. This is a direct consequence of the fact 
that each mechanism localizes the dynamics in a different 
sense and, when they are of the same order of the hop¬ 
ping amplitude, their competition partially cancels each 
other’s effect. This non-trivial interplay is also imprinted 
in the gap distribution of the Hamiltonian indicating that 
it also affects the integrability of the system. Finally, the 
interplay is shown to become more pronounced with in¬ 
creasing system size, robustly showing at least partial 
disorder-interaction cancellation. 
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Appendix A: Equilibration analysis 

Proposition: Let Ti. = (C^)®'^ be the state space of a iV sites quantum system with non-degenerate Hamiltonian 
H. Assume 1/ {{L)) < 1 -|- d , with (( )) = ^ where the sum runs over the elements |7r) of some fixed product 

basis H, L(|7r)) = J2e I (’’’ 1 -^) the sum running over all eigenstates \E) of H, and 0 < d < 1. For every 

subsystem S with Ns sites satisfying there exists a pair of orthogonal product states Itt^) and |7fs) on 

S and a product state Itta) on the remaining system R, such that 

\im ^ f D{ps{t),psit))dt>l-6V6, 

T->oo T Jq 

where ps{t) is the reduced state of S under the evolution of the system with initial state |7rs) G Ittr), ps{t) is the 
reduced state of S under the evolution of the system with initial state |7fs) G and D{ps{t),p{t)) = \\\ps — Ps||i 
is the trace distance. It is thus shown that if the dynamics of the global system is localized in state-space then it is 
guaranteed that large enough subsystems composed of Ns sites do not equilibrate. This means that correlations do 
not extend over Ns sites and thus the correlation length should be much smaller. 

Proof. Take a subsystem with Ns sites satisfying Assume, by contradiction, that for each product 

vector Itt^) on the remaining system R, there is at most one product state [tts) such that {its) ® Itt^) G H and 
L{\Trs) G ki^)) > I — <5. We would therefore have at least d^~^^ x (d^® — 1) elements of H satisfying L{\Tr)) < 1 — S. 


But then 


= = E ^(W)+ E ^(W)] 


7r:L(|7r))<l —(5 


7r:L(|7r))>l —(5 


< ^[(1 - 6).{d^^ - + l.d^-^^] 

<(l-<5)(l-y^)+ 


S + l' (5 + 1 


l + (5’ 


(Al) 

(A2) 

(A3) 

(A4) 

(A5) 


a contradiction. 

From the argument above we must have then two orthogonal product vectors [tt^) and |7rg) on the state space of S 
and a vector |7ri^) on the space of R such that L{\tts) 0 Ka)) >1 — 5 and L{\Trs) 0 |7r_R)) >1 — 5. 

Let |7r(t)) be the evolved state of an element [tt) of H. Since we can write, for every t > 0 


\Tr{t)) = a{t)\7T) + I3{t)\'ip{t)) and \n{t)) = a(t) \n) + P{t) , 

where a{t) = (7r|7r(t)), |■i/'(^)) is normalized and orthogonal to [tt) for all t (so |Q;(t)P + |/3(i)P = 1)) analogous 
dehnitions for |7r(t)), we have: 

D{ps{t),ps{t)) = D{Tyr |7r(t)) (7r(f)| ,Tri^ |#(t)) (if(<)|) (A6) 

= ^l|Tri?[|a(t)p Itt) (tt] + |/3(t)|2 \^p{t)) {^P{t)\ + ia{t)P{t)* K) {^P{t)\ + h.c.) (A7) 

- id(t)p 1^) (7f| - \m? \m) {m\- m)m iti) {m\ + h.c.)iii (as) 

= ^11 ks) ksl - ks) ks| - |/3(t)p ks) ks| + |/3(t)p ks) ksl (A9) 

+ Trfi[|/3(t)|2 \'ij;{t)) + {a{t)l3{t)* |7r) (V'(t)| + h.c.) (AlO) 

- \m^ m)) {m\ - {a{mtr k) m)\ + h.c.)]\\u (ah) 


were in the last equality we have taken the partial trace of projectors |7r) (7r| and Iff) (7f| and used that |a(t)p + |/3(t)P = 
1 and |d(t)P + |/3(i)P = 1- Using the reverse triangle inequality in the last expression, we have 


D{ps{t),ps{t)) > ^11 ks) (ti-sI - ks) (^sl 111 - ^11 - l/3(t)P ks) ksl + k(f)P ks) ksl 
+ Tri^[|/?(f)|2 (V'(t)l + |7r) (V:(t)| + h.c.) 

- \m^ m)) (k(t)i - {a{t)m* k) {m\ + h.c.)]\\Y 

> 1 - mt)f - - 2\ait)mt)\ - ^Hmm 

>i-mt)\-mt)\. 


Now, since the Hamiltonian is non-degenerate, we have the well-known relation: 


lim ^ j |a(t)|2 = lim ^ f \ (7r|7r(t)) p 

t^oo 1 Jq t^oo 1 Jq 

= lim ^ r \Y.\{7T\E)\^e-^^r 

t^OO 1 jp. ^^ 

= ^ r E I (^1^) I'l 

t->oo J /„ ^' 

= Ek^i^) 

E 

= ^(k)), 


E,E' 
|4 


(A12) 

(A13) 

(AM) 

(A15) 

(A16) 


(A17) 

(A18) 

(A19) 

(A20) 

(A21) 
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and similar one for L(|7f)). Moreover, by the Cauchy-Schwarz inequality we have 


lim — 
T^oo T 


|/3(t)|dt< lim - 

T—^oo 1 


J Idt 

= \/l - L{\tt)) < Vs, 


and similar inequality for /3(t). Finally, from Inequality (A16), we get 


1 

lim - / D{ps{t),ps{t)) > 1 - gVS. 
T-i-oo T Jq 


(A22) 

(A23) 

(A24) 
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